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We study the evolution of the optical conductivity in the t-J model with temperature and dop- 
ing using the Extended Dynamical Cluster Approximation. The cluster approach results in an 
optical mass which is doping independent near half filling. The transition to the superconducting 
state in the overdoped regime is characterized by a decrease in the hole kinetic energy, in contrast 
to the underdoped side where kinetic energy of holes increases upon superfluid condensation. In 
both regimes, the optical conductivity displays anomalous transfers of spectral weight over a broad 
frequency region. 



PACS numbers: 71.27.+a,71.30.+h 

How superconductivity emerges from a strongly corre- 
lated normal state is one of the most important problems 
in the field of strongly correlated electron systems. This 
problem has received intensive attention of the past two 
decades spurred by the discovery of high temperature su- 
perconductivity. A large number of experimental probes 
have been applied to this problem, but there is, at this 
point, no theoretical concensus on the basic physics or the 
minimal model required to describe this rich phenomena. 

In recent years, significant progress has been achieved 
through the development of the Dynamical Mean Field 
Theory (DMFT) Q and its generalizations HUH- This 
approach, allows us to isolate which aspects of the physics 
of a given material can be understood within a local ap- 
proach, and what is the minimal cluster size required 
to describe certain phenomena. For example it has been 
shown that a single site method describes well the physics 
of the finite temperature Mott transition in materials 
such as V2O3 in a broad region of temperature and pres- 
sures around the Mott endpoint, while the interplay of 
Peierls and Mott instabilities in materials such as VO2 
and Ti203 requires a two site link as a minimal refer- 
ence frame 0. There is an evidence that the Hubbard 
or the t-J model treated by a four site cluster DMFT ap- 
proach captures many essential aspects of the physics of 
the cuprates such as its phase diagram 0, 0, Q j the vari- 
ation of the spectral function and the electron lifetime 
along the Fermi surface 0,01, E3- This approach system- 
atizes and extends the early slave boson treatments of 
the cuprates which now appear as restricted low energy 
parameterizations in a more general approach. 

One of the most powerful bulk probes of carrier dy- 
namics is optical conductivity^ Previous studies of optics 
in t-J and Hubbard model 11, 12] have allowed an accu- 



rate description of many experimentally observed anoma- 
lies found in the cuprate superconductors at high tem- 
peratures, around or above the normal temperature. In 
this letter we use the cluster extension of DMFT, which 
allows us to study the evolution of the optical conduc- 
tivity in the most interesting temperature regime around 
the transition from anomalous metal to d-wave supercon- 
ductor. As a result, we gain several new insights and we 
show theoretically that: a) while optical spectral weight 
is reduced in the underdoped side, the optical mass of 
the carriers remains constant, b) The changes in hole ki- 
netic energy when entering the superconducting state are 
much smaller than in the Hubbard model, and are posi- 
tive in the overdoped side and negative in the underdoped 
regime, c) The broad redistribution of spectral weight 
which is observed in the superconducting state is largely 
due to the anomalous Greens function, and in this re- 
spect optical measurements combined with theory allows 
us to extract additional information not available from 
photoemission spectroscopy, d) At optimal doping, the 
low frequency part of optical conductivity is Drude like 
while the intermediate range shows approximate power- 
law w~ 2 / 3 . All these findings are consistent with exper- 
iments and can be understood in terms of a plaquette 
embedded in a correlated medium but not in terms of a 
single site theory. 

Formalism Our starting point is the t-J model which 
was proposed by P.W. Anderson as a minimal model to 
describe the cuprate superconductors. We use an exact 
reformulation of this model in terms of a spin fermion 
model using the Hubbard-Stratonovich transformation 
to decouple the spin interaction. The effective action de- 
scribing the interaction of fermions with spin fluctuations 
takes the following form 
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The Hubbard U term in the action will be taken to in- 
finity to enforce the constraint. 

The exact Baym-Kadanoff-like functional for this prob- 
lem is 

T[Q,V] = -Trlog(Go 1 - E) - Tr[£E] 

+ iTrlog^ 1 - n) + ^Tr[Pn] + $[0,2>](2) 

We chose Extended Dynamical Cluster Approximation 
method (EDCA) 0, 0, 113 on a plaquette because it is 
the simplest mean field theory that justifies ignoring the 
vertex corrections of optical conductivity. The only ap- 
proximation of the EDCA is to replace the Green's func- 
tions in the interacting part of the functional <&[G, T>] with 
the corresponding course grained cluster Green's func- 
tions £ k -> g K = J^keK Gk and ~* V Q = S qG Q V Q 
where the sum X^keK ^ s over those k momenta in Bril- 
louin zone which correspond to certain cluster momenta 
K (see 6]). The saddle point equations of the functional 
Eq. Hare E K = S$/8g K and Hq = 28<$>/8V Q which 
automatically give pice-vise constant self-energies. To- 
gether with the Dyson equations C?k = X)keK:( G o 1 ~ 
Sk) _1 and 2?q = 5Z qe q(-Dg 1 — Hq) -1 form a closed set 
of equations. Few comments are in order: i) The bosonic 
self energy is simply related to the spin susceptibility 
Xq = (^Hq 1 + Jq)^ 1 an d this expression can be 
used to calculate bosonic self-energy knowing local spin- 
susceptibility, ii) The Baym-Kadanoff functional $[C?,P] 
which depends only on the cluster Green's functions can 
be obtained by solving the cluster problem coupled to 
the fermionic and bosonic bath 0, El • As an impurity 
solver we use exact diagonalization for the cluster and an 
NCA approach to treat the hybridization of the cluster 
with the bath. To study superconductivity we allow the 
off-diagonal long range order by employing the Nambu 
method of anomalous Green's functions. 

The optical conductivity of the plaquette within 
EDCA is particularly simple because the vertex correc- 
tions vanish just like in single site DMFT. The reason is 
that all four cluster K-points of the plaquette are course- 
grained over the part of the Brillouin zone which is partic- 
ularly symmetric, namely k as well as — k point for each 
k point are in the same cluster region. Since the current 
vertex depends only on the cluster K point, the vertex 
legs can be closed and coarse-grained X)k£K w k{?k * G\t- 
Since velocity is even function of k and C? k at fixed cluster 
K depends of k only through ek, this quantity vanishes. 

The optical conductivity can therefore be simply ex- 
pressed by 
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Results - In Fig. ^ we show the evolution of the op- 
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FIG. 1: Optical conductivity across the superconducting 
transition for various doping levels from underdoped to over- 
doped regime. Letter S in legend stand for superconducting 
state and N for normal state of the t-J model. 



tical conductivity for various temperatures and dopings 
from room temperature T ~ 0.09i down to transition 
temperature T c w 0.034t and slightly below the transi- 
tion T w 0.76T C . We have representatives of the three 
regimes of the cuprate superconductors, slightly under- 
doped S = 0.14, optimally doped S = 0.16 — 0.2 and 
the overdoped case 5 = 0.26. In agreement with exper- 
iments on high Tc's, and with earlier theoretical single 
site DMFT studies, the optical conductivity extends over 
a broad frequency range, and consist of a broad Drude 
component at low frequencies that sharpens with decreas- 
ing temperature and a higher frequency (" mid infrared" ) 
component with substantial intensity at high frequencies. 

Upon entering the superconducting state no real opti- 
cal gap opens below the transition, unlike the standard 
case of the s-wave BCS superconductors. However, a 
substantial reduction of the optical conductivity due to 
superconductivity is observed up to very high frequency 
even beyond lo > t « 2500cm -1 which is ten times bigger 
that the superconducting gap at optimal doping A ss O.lt 
itself. 

The Drude-like low frequency conductivity narrows as 
the temperature decreases from room temperature to just 
above T c and continues to grow even below the transition 
temperature. This enhancement of the low frequency 
conductivity as a result of the onset of coherence is seen 
in many experimental studies and should be con- 
trasted with the reduction of the conductivity at higher 
frequencies due to the opening of a gap. 

The optical conductivity has a very long tail of incoher- 
ent spectral weight that dominates the optics. Further- 
more, the tail seems to have an approximate power law 
a cx (— iuj) a with the exponent close to a — 2/3 ^1 in 
agreement with experiments [l^. Hjj . At low frequency, 
optical conductivity is Lorentzian-like at optimal doping 
with w/T scaling. In Fig.[2t we plot (T<7(w)) -1 as a func- 
tion of u/T which is temperature independent quadratic 
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parabola in the region —T < 10 < T in agreement with 
recent experiments [l^ . 

The peak around 0.15i is the standard BCS coherence 
peak arising from the excitations across the supercon- 
ducting gap in one particle density of state and is not 
visible in experiment which is closer to the clean limit 
than our theoretical calculation which overstimates the 
scattering rate. 
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FIG. 2: a) Optical effective mass as determined from Eq. [I] 
and plasma frequency versus doping. The plasma frequency 
decreases towards the Mott point because of decrease of the 
active part of the Fermi surface while effective mass remains 
almost constant, b) The integral of the optical conductiv- 
ity J Q " cr'(x)dx for normal and superconducting state. The 
optical weight in SC state missing in the finite frequency re- 
gion collapses into 5 function. The spectral weight is restored 
when integral is taken up to ~ At ~ leV ~ 8000cm -1 . c) The 
inverse of low frequency optical conductivity l/(Ta(ui)) ver- 
sus cu/T is quadratic parabola independent of temperature 
in the optimal doping regime. Optical conductivity is thus 
Drude-like at low frequency up to ui ~ T. 

Optical conductivity measurements can be used to de- 
fine an optical effective mass defined by 
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The plasma frequency in this equation is proportional to 
the kinetic energy of the model and at the same time, to 
the integral of the optical conductivity. This definition 
can be motivated by the extended Drude model of ref- 
erence^^. Here we see that the cluster DMFT method 
introduces new physics not captured by the single site 
approach. 

In single site DMFT the optical mass diverges as the 
Mott transition is approached. Taking into account the 
low frequency and low temperature expansion of the self- 
energy in the local Fermi-liquid state = fJ, — Ho + 
u>(l — ^) — iT one obtains for the conductivity a(ui) = 
e 2 &xx([J<o)(l + iuj/(2ZT))/(2T) and therefore the optical 



mass in single site DMFT is m*/m = 1/Z and diverges 
as Mott transition is approached since Z vanishes. 

In the cluster case, the mass no longer diverges. In 
Fig. we show the optics mass from Eq. 0] where oj p 
was computed by integrating the optical conductivity up 
to At w leV. When approaching Mott transition, the 
effective mass becomes 5-times bigger than the electron 
mass while it drops down to 3-times the electron mass in 
the overdoped regime. This behavior is very similar to 
the recent experiments on cuprates reported inl20t The 
plasma frequency in Fig. which is proportional to the 
integral of optical conductivity, grows linearly at small 
doping and the slope gets less steep around 10% doping 
in accordance with recent experiments on cuprates |2lj . 

The EDCA equations in the 2x2 plaquette provide 
us with a simple description of how the the Mott insu- 
lating state is approached, since one has only four differ- 
ent patches in momentum space, which however, change 
very differently with doping. The self energy that is 
relevant to the optical conductivity corresponds to the 
patch around (0, 7r) and (tt, 0) and has a quasiparticle 
residue which is approximately doping independent and 
finite (Z k 0.1) as the Mott insulator is approached. 
However, the optical conductivity vanishes with vanish- 
ing doping just like the plasma frequency in Fig. [3 This 
effect is driven by the vanishing of the number of car- 
riers measured by the length of the active segment of 
the Fermi surface, which goes to zero as the Mott point 
is approached, rather than with vanishing of quasiparti- 
cle renormalization amplitude. This interpretation of the 
EDCA results is consistent with the results obtained by 
the cumulant periodization of CDMFT [22j |. 

Transfer of spectral weight - A surprising aspect of the 
physics of strongly correlated materials, is that low en- 
ergy phenomena affect the spectra of the material over 
a very large energy scale. This general phenomena is il- 
lustrated in Fig. |2Jd, which shows the integral of optical 
spectral weight in the normal and the superconducting 
state. Low energy phenomena like the onset of supercon- 
ductivity which involves a scale of a fraction of J, involves 
redistribution of optical weight of the order of At ~ leT^ 
which is 40 times more than the gap value. This large 
range of redistribution of spectral weight was also mea- 
sured on cuprates and pointed out in |U A new 
theoretical insight is that the high frequency redistribu- 
tion of weight comes from the anomalous Greens function 
T * T in Eq. @ and hence can not be observed in the 
density of states or ARPES measurements. 

Sum Rule - To understand the mechanism of high tem- 
perature superconductivity, much effort was recently de- 
voted to measure the chan ge o f kinetic energy upon su- 
pcrfluid condensation (2J, |2jj|. Experimentally this is 
measured by integrating conductivity up to large enough 
cutoff of the order of leV in both normal and supercon- 
ducting state, and interpreting the results in terms of the 
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optical sum rule which states 

f A ire 2 
J o <j'{x)dx = - — {T) (5) 

This sum rule has also been extensively discussed in the 
context of the Hubbard model, in which case, assum- 
ing nearest neighbor hopping matrix elements only, the 
operator on the right hand side of Eq. [SJis the kinetic en- 
ergy of the electrons in the Hubbard model. It has been 
shown by Jarrell and collaborators [2|| that entering the 
superconducting state results in a reduction of kinetic 
energy, in stark contrast with the BCS model where the 
kinetic energy increases upon entering the superconduct- 
ing state. 

The kinetic energy of the Hubbard model is composed 
of two different contributions. The superexchange en- 
ergy of the spins, and the kinetic energy of the holes. 
In the Hubbard model the optical transitions at ener- 
gies of order U, give rise to the superexchange energy of 
the spins, while the low energy transitions correspond to 
the kinetic energy of the holes, but it is not possible to 
separate those two contributions in the Hubbard model. 
Here we sharpen the analysis of Ref. |2(| by evaluating 
separately these two physically different contributions, 
which allows us to make direct contact with optical ex- 
periments. Namely, the experiments measure the kinetic 
energy of the holes since the cuttoffs which are used are 
such that the transitions into the upper Hubbard band 
are not included. 




0.028 0.03 0.032 0.034 0.028 0.03 0.032 



T/l T/l 

FIG. 3: The difference between the superconducting and nor- 
mal state energies as a function of temperature. The following 
curves are shown: blue with triangles up - E^in-s — Ekin-N', 
red squares - E exch _ s - E exch _ N ; green diamonds - E to t-s - 
Etot-N; black with circles - [is — fJ-N- 



we see that the change in kinetic energy upon superfluid 
condensation changes sign between overdoped to under- 
doped regime. In overdoped regime, the conventional 
BCS picture is applicable where kinetic energy of holes 
increases while the superexchange energy decreases just 
like the interaction energy in conventional phonon medi- 
ated superconductors. The later change is much larger 
such that the total condensation energy is positive. In 
underdoped regime, the superexchange energy still de- 
creases upon condensation and gives the largest contribu- 
tion to the condensation energy Q. However, the kinetic 
energy of holes is gained in uderdoped case so that the 
kinetic energy gives positive contributes to condensation 
in agreement with recent experiment |25|. 

In Fig ©, we also show the change of chemical po- 
tential upon condensation. Similarly to kinetic energy, it 
also changes sign with doping and it increases in under- 
doped regime while it decreases in overdoped regime. We 
believe this is a consequence of the asymmetry of the nor- 
mal as well as the superconducting one particle density 
of states. 

Conclusion - In conclusion, we have studied the t-J 
model using the EDCA in a plaquette near the optimally 
doped regime. The theoretical results capture the main 
features of the optical conductivity near the transition 
from a strongly correlated metal to a superconductor. 
In spite of the shortcomings of the approximation used 
(the EDCA 2x2 cluster overestimates the superconduct- 
ing critical temperature and the quasiparticle scattering 
rates) we believe that the lessons from this mean field 
theory, as to the physical origin of the scaling of the op- 
tical mass with doping, the anomalous redistribution of 
spectral weight and the changes of kinetic and exchange 
energy upon entering the superconducting state are gen- 
eral, and can be derived with considerable more effort 
in other cluster schemes, or using larger clusters. They 
stress the importance of the differentiation of states in 
momentum space in strongly correlated materials near 
the Mott transition. 
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In Fig. |3 we show the change of kinetic energy of 
holes, the superexchange energy (J2ij Jij '(SjSj)/2) and 
total energy upon condensation as obtained from our 
EDCA cluster calculation. The kinetic energy can be 
directly calculated from the spectral function and the 
superexchange energy from the spin susceptibility. In 
agreement with very recent experiments on cuprates |25| 
and in agreement with RVB and slave-boson prediction, 
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